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RELATIONSHIPS  AND  TRANSFORMATIONS 
BETWEEN  CONCENTRATION-PATH-LENGTH  (CL), 

AGENTS  CONTAINING  PARTICLES  PER  LITER  OF  AIR  (ACPLA),  AND 
THE  NUMBER  OF  SPORES  (  N^pom  ) 


1 .  INTRODUCTION 

In  remote  sensing,  passive  Infrared  (IR)  sensors  (e.g.,  Fourier-transform, 
FTIR  spectrometers),  measure  mass-column-density  (Concentration-path-Length  (CL) 
mg/m  ).  Aerosol  mass  concentration  is  represented  by  C  (mg/m  )  and  L  (m)  is  the  path- 
Length  along  the  line-of-sight  of  the  sensor  containing  the  aerosols.  The  threat  level 
(detection  objective)  is  usually  alluded  to  as  Agent-Containing-Particles  per  Liter  of  Air 
(ACPLA).  Agents  can  be  either  spores  (microbes,  toxins,  viruses)  that  are  “stuck”  together 
to  form  an  aggregate  ACPLA  particle  or  the  spores  can  be  cohered  to  the  surface  of  a  large 
particle  such  as  dust.  The  quantities  of  interest  are  listed  below: 

(a)  The  CL  (the  cumulative  mass  of  the  spores  along  a  detection  path-Length). 

(b)  The  number  of  spores  (N  es)  along  a  path-Length. 

(c)  The  ACPLA  (the  number  of  aggregates  in  a  liter  of  air). 

The  3  concepts,  CL ,  N spom ,  and  ACPLA,  are  of  interest  to  different  members 

of  the  community  for  different  reasons.  For  the  engineer  who  builds  an  IR  sensor,  CL  and, 
therefore,  N  s  and  their  cumulative  mass,  are  crucial  parameters  for  determining  the 

sensitivity  of  the  sensor.  The  relationship  between  N  es  and  CL  is  also  useful  for 

computing  fluorescent  cross-sections  from  bioaerosols.  For  the  person  who  designs  and 
evaluates  the  efficiency  of  delivery  systems,  ACPLA  is  of  interest.  A  sensor  based  on 
counting  particles  by  light  scattered  from  a  particle  also  uses  ACPLA.  For  a  physician,  the 
concept  of  N  s  is  critical. 


This  study  explores  the  relationship  between  CL ,  N spores ,  and  ACPLA ,  and 

illustrates  a  computational  procedure  that  transforms  ACPLA  to  CL  andAvOTra.  The 

difference  between  ACPLA  and  CL  is  determined  by  the  spore  count  in  an  aggregate  of  a 
certain  size,  which  is  dependent  on  bacteria  type  and  the  dissemination  scenario  for  the 
spores,  which,  in  turn,  gives  rise  to  the  issue  of  packing  density. 


2.  NUMBER  AND  SIZE  VOLUME  DISTRIBUTION 

The  aerosol  number  size  distribution  density  function,  f„(r) ,  is  modeled  as 

a  lognormal  distribution.  In  the  context  of  this  study,  fn(r)  will  be  regarded  as  the  size 

distribution  of  the  aggregate  of  size  (  r ).  For  convenience,  ( r )  will  be  regarded  as  the 
“radius”  of  a  spherical  aggregate.  The  aerosol  number  size  distribution  density  is  given  as 
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(1) 


/„(r  ://,cr)  = 


(lnr-ln//)2 
ex p[ - — — 77H 


^2^(lncr)2  '  ‘  2(lncr)‘ 

w  =  f/>  :  //,  cr)t/r  ( number  of  aerosol  particles  within  a  unit  volume ) 


where  p  and  o  are  the  parameters  of  the  log-normal  size  distribution;  f„(r)dr  is  the 

dr 

number  of  aerosols  within  radii  r  ±  — ;  and  N  is  the  total  number  of  aerosol  particles  in  a 

unit  of  volume.  When  the  unit  of  volume  is  taken  as  I  liter ,  ACPLA  =  N  ( liter  ' ) . 

The  equation  for  volume  size  distribution  density  function,  /v(r) ,  is  given  as 


fv(r:p,a)  =  fn(r\p,a) 


4m' 


(2) 


dr 

where  ffr)dr  is  the  volume  of  aerosol  particles  of  radius  r  ±  — .  The  total  volume  is 
V (p r,cr) .  Equation  2  and  all  aerosol  radii  sizes  are  integrated  to  determine  V  and  the  CL  : 


EO,cr)=  \^^fn(r)dr=  J 


4/zrJ 


r^/27r(ln  <j)‘ 


rexp[- 


[In  r  -  In  pf 
2(ln<7) 


T-\dr  =  \ffr)dr 


(3) 


3.  PARAMETERS  FOR  NUMBER  DENSITY  ffr)  AND  VOLUME 

DENSITY  ffr)  DISTRIBUTIONS 

The  parameters  that  describe  the  lognormal  distribution  f(r  :  p,cr)  are  mode 
radius  ( rmode ),  median  radius  ( rmedmn  ),  mean  radius  (  rmean  ),  and  standard  deviation  ( std  of 
the  distribution).  All  4  parameters  are  a  function  of  2  variables  (p,cr)  that  characterize  the 

lognormal  distribution  given  in  eq  1.  The  rmode  is  used  when  size  distribution  is  at  its 

00 

peak:  f(r  =  rmode)  =  max  .  For  rmedjan ,  defined  as  jf(r)dr  -  0.5 , 50%  of  the  distribution 

r  median 

CO 

area  is  above  the  radius.  The  rmean  parameter  is  defined  as  rmsaa  =  \rf{r)dr  (the  first 

o 

moment  of  the  distribution).  The  std  of  the  distribution  is  defined 

as  std  =  j(r  -  rmean )2  f(r)dr  with  a  lognormal  distribution,  std  *  a  .  While  rmnw , 

r,nedianand  std  are  all  functions  of  the  integration  limits  [rmjn,  rma  J ,  which  are  implicitly  0  and 
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co  in  eqs  4  and  5,  rmode  does  not  depend  on  the  radii  range  [rmin  ,rmax  ]  used  in  the  size 
distribution  function, /(r) . 

For  the  number  density  size  distribution  fn(r  :  /v,cr),  where  the  radii  (r)  is 
unbounded  (0  to  oo),  the  parameters  are 

rmea„  =exp(ln/r  +  0.5(lna)2) 

^mod  e  =  exp(ln/y  -(lna)2) 

<  > 

^  median  ft 

std  -  -y/exp[21n /7+  (In  cr)2  ][exp((lncr)2)  -  l] 

For  the  volume  density  function  fv(r  :  /j,  a) ,  the  parameters  are 


rmean  =exp[/  +  0.5(lno-)2] 
rmod,  =exp[/-(lncr)2] 

^ median  =  exp(/) 

std  =  yjexp[2t  +  (ln<r)2  |[exp((lncr)2)-  l] 
where  t  =  In  fj.  +  3(ln  cr)2 


(5) 


mode  mode 


Figure  la.  Number  and  Volume  Density  Lognormal  Distributions  for  a  Typical  Bioaerosol 

{H  =  0.875  /im,  a  —  2.09  uni). 

In  Figure  la,  number  and  volume  density  lognormal  distributions  are  shown 
for  typical  bioaerosol  parameters.  For  the  number  density  distribution,  f„(r),  the 
parameters  are 


W  =  0-5^ » 

Median  =  0-87/7/77  , 
rn,ean  =  1  - 1  4/7/W  ,  and 

std  =  0.91  fum  . 

For  the  volume  density  distribution, /,(/) ,  the  parameters  are 
rmod,  =  2.58/vm , 

^median  =  4  44 V™  > 
rmean  =  5-83/7/77  ,  and 

std  =  4.95 jum . 

For  both  Figures  la  and  lb,  the  parameters  for  the rmcdlan ,  std,  and  rmcun , 
which  are  dependent  on  the  value  of rmax ,  are  for  size  distribution  with  radii 
0  <  r  <  (rmax  — »  oo) .  The  location  for  rmode  is  independent  of  rmax . 


Figure  lb.  Number  Density  f„(r:  //,  o)  and  Volume  Density  fv(r:  /u,  a)  Lognormal 

Distributions  (//  =  1  ,um,  o  =  2  p.m). 

The  lognormal  distributions  for  number  density,  fn  (r  :  //,  a) ,  and  volume 
density,  fv(r  :  ,u,a) ,  are  illustrated  in  Figure  lb.  For  convenience,  this  size  distribution 
was  used  to  represent  bioaerosol  throughout  the  examples  given  in  this  study  (Figures  la 
and  lb  are  very  similar).  The  parameters  for  the  number  density  distribution,  f„(r) ,  are 
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'mod,  =  °-62 ,Lim  , 

Median  =  >/"«  , 

rmean  =  1.27/»,  and 
std  —  0.99 /jm  . 

The  parameters  for  the  volume  density  distribution  /v(r)  are 

'mode  =  2.61pm  , 

'■»«**  =  4.23/rm , 

'■*«»  =  5.37pm ,  and 
s/r/  =  4.22 pm  . 

Due  to  gravitational  settling,  some  particles  fall  to  the  ground  while  traveling 
a  distance  (x)  with  a  wind  velocity  (w)  within  the  time  frame  (t=x/u).  The  falling  particles 
can  also  be  subjected  to  other  deposition  processes  such  as  turbulence  to  a  non-reflecting 
surface  and  rain.1  When  a  particle  of  radius  (r)  is  dropped  from  a  height  (//),  the  falling 
velocity  will  at  first  increase  (acceleration  by  gravity)  and  then,  after  a  time  (t)  and  a 
distance  ( h )  (due  to  the  drag  force),  the  particle  will  reach  a  terminal  settling  velocity  for  the 
remaining  distance  to  the  ground. 


0  1 


Figure  2.  Terminal  Settling  Velocity  of  Sphere. 


Figure  2  portrays  the  terminal  velocity  of  a  particle  of  radius  (r)  and  density' 
of p  =  1.32  gem ~3  falling  in  still  air  at  20  °C  temperature.  The  time  taken  for  reaching 
terminal  velocity  is  less  than  0.01  s.  A  particle  with  a  radius  of  25  /urn  ,  dropped  from  a 
height  of  100  m  and  carried  out  by  a  wind  at  a  velocity  of  3  m/s  (~10  km/hr),  will  travel  a 
distance  of  3  km  before  it  hits  the  ground.  Thus,  great  care  is  needed  in  choosing  the 
maximum  radius  in  a  given  aerosol  size  distribution  for  a  particular  scenario. 
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4. 


TOTAL  VOLUME  AND  UPPER-LIMIT  CL  ESTIMATE  FOR  AEROSOL 
SIZE  DISTRIBUTION 


-r  10 


01  23456709  10 

M  0‘<n) 


Figure  3.  Total  Volume  Vfa,  a)  of  All  Particles  in  Lognormal  Number  Density  Distribution 

f„  (r:  jj,  a). 


In  Figure  3,  the  total  volume,  V{p,a),  as  discussed  in  eq  3,  of  all  the 
particles  is  shown  in  the  lognormal  number  density  distribution,  fn{r  :  The  total 

number  of  particles  (N),  in  fn(r) ,  within  a  unit  volume  is  assumed  to  be 

^max 

N  =  J f„(r)dr  =  1 .  The  range  of  particle  radii  is  between  rmin  — >  0  and  rmax  =  25  fjm . 


When  the  unit  of  volume  is  taken  to  be  1  liter,  ACPLA  =  N  ( liter  1 )  . 

10* 


Figure  4.  Same  as  Figure  3,  but  the  Largest  Radius  ( rmax  =  5/jm )  Shows  Only  Small 

Particles  Being  of  Interest. 
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Figure  4  also  shows  total  volume  V(p,a) .  The  largest  radius  (rmax  =  5pm ) 

reflects  a  scenario  in  which  only  small  particles  are  of  interest.  Some  of  the  curves  are 
shown  to  cross  each  other  around  p  ~  2.5  pm  .  This  seemingly  strange  behavior  is  ascribed 

to  the  fact  that  for  each  pair  of  ( p,G ) ,  the  total  number  of  particles  (N)  in  fn(r)  is  the  same. 
The  median  and  mode  of  the  distribution  is  a  function  of(p,a) .  As  the  value  of  (//,cr) 
increases,  most  particles  (mode  of  distribution)  shift  toward  a  larger  radius:  for  p  ~  2.5  pm 
and  g  >  1 .15/jm  ,  rmode  shifts  to  r  >  5  pm .  However,  due  to  the  constraint  of  rmax  =  5 urn  in 
the  size  distribution,  particles  with  r  >  5pm  (the  expected  mode  of  distribution)  cannot  be 
attained.  Thus,  the  increase  in  overall  volume  is  possible  only  by  the  small  particles  with 
r  <  5 pm  whose  relative  number  is  a  function  of  (/ j,g ) .  The  rate  of  increase  for  curves 
with  smaller  a  is  faster  due  to  the  fact  that  the  mode  of  distribution  has  not  yet  been 
reached. 


The  conversion  of  total  volume  (V)  (the  y-axis  of  Figures  3  and  4)  to  CL  is 
straightforward.  The  unit  volume  assumed  for  the  total  number  of  particles  (N)  in  fn(r) 

can  be  taken  to  be  1  liter  (1000  cm  ).  Thus,  for  a  given  aerosol  density,/?  (— =— )  and  a 

cm 

path-Length  of  I  (m) ,  the  Concentration-path-Length  product  (mass-column-density,  CL)  is 


,mg  particles  5  liters  3 
CL(—=r)  =  N- — : - xio  — —  V (cm  )x/>(- 


m 


liter 


m 


cm 


r)x  103  — —  x  L(m)  (6) 

g 


The  CL  in  eq  6  symbolizes  the  aerosol  material  (/?)  in  the  size 
distribution  fn{r  :  ju,g).  Not  all  the  aerosol  mass  in  fn{r  :  p,o)  consists  of  the  desired 

agent  material  ( BG  spores)  because  of  the  packing  efficiency  (the  fractional  volume  of  the 
spores  within  the  aerosols)  of  the  spores  into  aggregate  “particles”.  Thus,  the  computed  CL 
value  in  eq  6  may  be  regarded  as  an  upper-limit  estimate  whereas  the  correct  CL  will  be 
much  lower  when  the  volume  occupied  by  the  spores  in  the  aggregate  is  accounted  for 
(Section  8). 

The  N  s  present  in  the  given  size  distribution,  f„(r) ,  is  obtained  by 
calculating  for  CL.  For  1000  ACPLA,  a  100-meter  path-Length  and  a  size  distribution 
fn{r  .  p  =  \,g  =  2)  (Fig.  lb)  with  (BG)  a  density  of p  -  1.32  gl cm3 ,  the  (upper-limit) 
estimate  for  CL  derived  with  eq  6  is 

CL  =  1000( )  X  4  X  1 0  1 1  (e/77 3 )  X  1 . 3 2  (g/cm3)  x  100(/n)x  106  =5.3 

liter  m~ 


for  particles  with  r  <  25 pm  (Figure  3).  When  only  small  particles  with  r  <  5pm  (Figure  4) 
are  considered,  CL  will  be  calculated  as 
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CL  =  1 000 ( Particle) x  2 x  1 0  "  (cm3 )  x  1 .32  (gW)x  100(m)x  106  =2.6  ~ . 

liter  m ' 

Density  of  Bacillus  subtilis  BG  spores. 

Only  a  few  references  to  the  density  of  Bacillus  subtilis  BG  spores  are 
available  for  light  and  heavy  spores:  1 .45  g/cm \2  1 .2  g/an  ,3  1 .290  g/cm3,  and  1 .335 
g/cm  .  '  The  difference  between  the  values  may  partly  be  a  result  of  the  amount  of  water  in 
a  BG  particle.  Apparently,  no  “hard  numbers,”  which  are  extremely  important  for 
estimating  mass-column-density  CL,  exist  for  BG  density.  In  this  study,  the  average  of  all 
the  reported  densities  was  taken  to  attain  the  value  of  1 .32  g/cm  . 

5.  PHYSICAL  DESCRIPTION  AND  PROPERTIES  OF  BACTERIA  AND 

SPORES 

Typical  members  of  several  bacterial  species  are  shown  in  Figure  5. 5  Each 
of  the  first  top  3  bacteria  is  represented  by  a  single  spore.  The  last  3  bacteria,  the  bacterium 
being  the  spore,  are  taken  to  be  single-celled  organisms.  The  bacteria,  modeled  in  a  super¬ 
ellipsoid  asx-^-z  coordinates,  are  given  as 

(£)"+(T)" +(-)”=  1  (?) 

a  b  c 

The  semi-axes  of  the  ellipsoid  are  a ,  b  and  c.  When  a  =  b  <  c  and  n  =  2,  eq  7  describes  the 
regular  prolate  and  when  a  —  b  =c,  n  =2,  eq  7  is  a  sphere  with  radius  a. 
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Bacillus  anthracis 


Bacillus  atrophaeus 


| _ |  Vegetative  Cell 

■  Endospore 


Bacillus  anthracis  (SPECIAL) 


O 

Staphylococcus  aureus 


O 

Staphylococcus  epidermidis 


m 


Escherichia  coli 


o  1 

MICRONS 


Figure  5.  Typical  Members  of  Several  Bacterial  Species5  where,  for  the  First  Top  3 
Bacteria,  a  Single  Spore  is  Inserted.  The  last  3  bacteria  are  taken  to  be  single-celled 
organisms  (the  bacterium  is  the  spore). 

The  volume  of  the  super-spheroid  is  calculated  as 


V 


8  a  be 

? 

n~ 


/?(- +i,  -)/?(- ,-) 
n  n  n  n 


(8) 


l 

The  mathematical  beta  function  is  f3(z,w)  =  fr'”1  (1  -  t)w~'  dt  and  the  volume  equation  6  has 

o 

been  modified  for  eq  7. 5 

Table  1  contains  physical  parameters  for  the  bacteria  and  spores.  Additional 
reports  on  elliptically  shaped  (~1  long  and  ~0.7  ^  wide)5  Bacillus  subtil  is  BG 
spores7"9  and  rod-shaped  (-0.5  ^ m  wide  and  1-2^™  long)  Erwinia  herbicola]0  ] 2  are  in 
existence.  Estimating  the  Nsporev  is  crucial  for  estimating  the  CL  that  can  be  contained  in  a 

particle  of  radius  (r).  This  issue  is  addressed  in  2  different  ways  in  Sections  6  and  8. 
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Table  1.  Physical  Properties  of  Selected  Bacteria  and  Spores. 
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6. 


Nspons  AND  CL  ESTIMATE  VIA  A  “RULE  OF  THUMB 


A  widely  used  rule  of  thumb  propounds  that  the  number  of  spores  contained 
in  a  particle  (aggregate)  of  radius  (r  in  microns)  is  (2r)2  55 .  This  rule  also  states  that  a 
particle  with  r  <0.5  pm  cannot  contain  a  spore  and  a  particle  with  r  =  0.5  ptn  is  a  single 
organism  (spore).  For  the  lognormal  number  size  distribution  fn(r  :  p,c),  the  total  TV  t 
within  a  unit  volume  and  the  CL  are  calculated  as 

1  ma* 

Nspores  =  N  j  (2r)2-5i  /„  (r  :  p,  cr )dr 

0 .5 /.an 

■CL=(N,rm,^)x(p,m,V^r,)xL 

N  =  ACPLA  =  T  f„(r :  fi,(r)dr  (Partlcie*  -) 

J  liter 

rmm 

where  pspon,  is  the  density  and  Vspore  is  the  volume  of  a  single  spore  (eq  8)  so  that 
PsporeV sport  =  m spore  's  the  mass  (ms)  a  single  spore.  The  Nspores  x  10  !  are  the  number  of 

spores  within  a  cubic  meter,  ( sPores ) ;  an(j  i  js  the  path-Length  (m).  The  integration  in 

m 

eq  9  is  started  at  r  -  0.5 pm  ,  which  is  the  particle  with  the  smallest  size  radius  within  the 
framework  of  the  rule  that  can  contain  a  single  spore. 


spores 

liter 


<») 


(9) 


0  05  1  1  5  2  25  3  35  4  45  5 

(j-im) 

Figure  6.  Nspores  Using  “Rule  of  Thumb”  Eq  9  in  the  Lognormal  Number  Density 
Distribution  fn{r  :  p,<r)  with  Particles  =  r  >  5 pm . 
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In  Figure  6,  the  Nsporcs  estimate  for  a  family  of  lognormal  number  density 
size  distribution  of fn(r  :  p,cr)  is  shown  with  ACPLA  =  1  liter .  The  radii  range  in 
fn(r  :  ju,a)  is  rmm  ->  0  and  rmax  =  25  pm.  While  rmin  ->  0  in  Figure  6,Nspom  is  computed 
only  for  r  >  0.5 pm  (eq  9).  Figure  13  shows  that  for  the  lognormal  size  distribution 
/n(r  :  //  =  l,cr  =  2)  (Figure  lb),  the  desired  value  for  a  1-  A/'  f  particle  is  20  (Ar=  1 

particle  in  Figure  13).  The  estimated  CL  for  BG  spores  for  1000  ACPLA ,  and  100  meters 
path-Length  is  given  in  eq  9  as 


spores 

cl= mCarliclesx20  ^ 

liter  particles 
liter 


xlO^xO.356-10-9 

m 


mgx  1 00m  =  0.7 1 

m 


The  particle  size  distribution  fn{r.p  =  \,o  =  2)  contains  only  particles  smaller  than5/v/?7 
and  =0.356-10-|!g^l0j^-  =  0.356-10-*mg  (Table  1). 


The  mass  of  the  £G  spores  (0.356  10  l2g),  not  the  mass  of  the  BG  bacterium 
(2.93  10"12  g ),  was  chosen  on  the  assumption  that  only  the  spore  material  contains  spectral 
information  for  detection  algorithms.  The  bacteria  material  is  mainly  water.  Assuming 
both  bacteria  and  spore  density  are  the  same  (1.32  g/cm  }),  the  CL  from  the  bacteria,  l  spore 
per  bacterium  (20  bacteria/1  for  liter  particl e/liter),  is  much  larger  due  to  mass  ratio: 


0.712x 


2.93 

0.356 


=  5.86^. 


m 


7.  PACKING  EFFICIENCY  FOR  BACTERIA  AND  SPORES 


Packing  efficiency  of  spheres  (the  fraction  of  volume  filled  by  spheres)  is  an 
old  unresolved  issue.  In  1611,  Kepler  postulated,  but  did  not  prove,  that  the  maximum 

n 


packing  efficiency  of  spheres  is  ^ 


Kepler 


3V2 


=  0.74  (“Kepler  problem”).  Over  350  years 


later,  the  theory  has  finally  been  proved  by  N.J.A  Sloan.14  The  packing  of  spheres,  cubes, 
ellipsoids  and  cylinders  in  a  wide  array  of  lattices  (cubical,  hexagonal,  tetrahedral) 
continues  to  be  a  current  research  topic  and  many  papers  are  routinely  published  on  this 
subject.  The  experimental  and  theoretical  aspects  of  packing  have  been  repeatedly 
addressed  by  mathematicians  and  businesses  (packing  medical  capsules  and  grains). 


The  method  of  packing  particles  is  directly  related  to  size  distribution  and 
shape.  Particles  can  be  packed  in  many  ways:  random  packing,  pouring  particles  one  by 
one  under  gravity  till  the  system  becomes  congested,  shaking  the  container,  etc.  The  list 
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seems  never  ending;  however,  this  study  deals  specifically  with  spheroid  and  rod  shaped 
spores  and  particles. 


length 

The  aspect  ratio  77  =  — - —  of  particles  (Table  1)  and  size  variation  play  a 
width 

major  role  in  packing:  “a  short  spherocylinder  (rod)  that  may  not  fit  in  an  interstice  when 
oriented  in  a  given  direction  may  fit  when  the  orientation  is  changed”.1 5  Packing  efficiency 
can  be  studied  through  the  interplay  between  the  excluded  volume  produced  by  caging  a 
particle  in  contact  with  other  particles  and  the  number  of  contacts  required  to  cage  a 
particle.  Williams  and  Philipse  qualitatively  explain  packing  efficiency:  “for  [a]  very  short 
spherocylinder,  the  excluded  volume  depends  only  weakly  on  the  aspect  ratio  as  opposed  to 
the  number  of  contacts  necessary  to  cage  the  spherocylinder.  As  the  aspect  ratio  increases, 
more  contacts  are  required  to  cage  the  particle  requiring  a  higher  volume  fraction.  [At] 
higher  aspect  ratios,  the  required  number  of  contacts  plateau;  however,  the  excluded  volume 
increases  significantly  and  drives  the  volume  fraction  (packing  efficiency)  back  down”.15 
The  packing  efficiency  of  a  long  thin  rod  is  inversely  proportionate  to  the  aspect  ratio  and  s 

given  as  <f>  =  — —  for  77  >  15  . 

2/7 

Table  1  illustrates  that  the  aspect  ratio  for  this  study  is  small,  77  <  15  .  Thus., 
the  above  formula15  seems  redundant  and  a  method  developed  specifically  for  the  packing 
of  spheroids  is  used.16  For  rod-shaped  bacteria  such  as  Erwinia  herbicola  and  Escherichia 
coli,  an  alternate  method  has  been  proposed.  The  model  fit  in  eq  10  for  the  data  on  packing 
spheroids16  is  in  reasonable  agreement  to  that  presented  in  Evans  and  Ferrar  for  rod  shaped 
bacteria. 


Figure  7.  Packing  Efficiency  as  a  Function  of  the  Aspect  Ratio  rj  =  —  of  Spheroids 

a 


(model  fit  calculated  with  eq  10). 


16 
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In  Figure  7,  packing  efficiency  as  a  function  of  the  aspect  ratio  for  packing 

x  o  y  o  z  o 

spheroids  is  (—)'+(—)“+(—)'  =1  (eq  7)  with  a  small  aspect  ratio,  <15.  Maximum 
a  a  c 

packing  efficiency  occurs  at  7  =  1.4  for  a  prolate  spheroid  and  r/  =  0.7  for  an  oblate 
spheroid.  The  approximate  symmetry  in  the  packing  efficiency  between  prolates  and 
oblates  <f>{rf)  =  ^(7"' )  is  clearly  shown.  A  model  (polynomial  of  third  order)  was  fit  to  the 
packing  efficiency  data  given  as 


<P  = 


0.3849  +  0.0893  Infa)  -  0.0999  ln(7)2  +0.0161  ln(//)3 
0.3849  +  0.0893  ln(|//|)  -0.0999  ln(|/7|)2  +  0.0161  ln^)3 


1  <  7  <  15 

1  . 

—  <  7  <  1 

15 


where  7  =  — 
a 


(10) 


The  values  of  the  packing  efficiency  tf>  for  bacteria  and  spores  are  listed  in  Table  1 . 


8.  CL  AND  Nspores  ESTIMATE  VIA  PACKING  EFFICIENCY 

The  value  of  the  packing  efficiency,^,  depends  on  the  aspect  ratio,  7,  which 
is  different  for  a  spore  (single  organism)  and  a  bacterium  containing  one  spore.  Therefore, 
2  scenarios  are  addressed  for  the  dissemination  of  biomaterial: 

Scenario  /:  bacteria  (vegetative  cells)  are  disseminated  and  consequently  aggregated  to 
form  a  size  distribution  of  fn (r  :  //, er) .  This  scenario  is  applicable  to  non-sporulating 
material  such  as  E.  coli.  Of  course,  the  number  of  bacteria  in  the  aggregate  equals  the 
number  of  spores.  The  y-axis  is  labeled  “  N  s.” 

Scenario  IE  only  spores  are  disseminated  and  consequently  they  aggregate  to  form  a  size 
distribution  of  fn  (r  :  /u,cr) .  This  scenario  is  applicable  to  sporulating  material  such  as 
Anthrax. 


The  N  s  and  CL  for  the  2  scenarios  are  given  as 
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V40lb**ri*m) 
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hfU'I&t  hi 


scenario  I 


N,r""-  1  >m„) 


v 

spore 


bacterium 

spore  >  it 

scenario  II 


spores 
liter 


,  mg 


volume  of  particles 


CL  =  (N  10  ) x(p  V  )xL  (— t) 

v  spores  '  \r  spore  spore  /  V  2  ^ 

m 

V{p,  a)  =  J  (r  :  //,  cr)Jr  ( 

J  3  liter 

rmm 

V,ACriA  =  ']f,(r:^Vr(l^) 

J  liter 


) 


(hi 


where  pspore  is  the  spore  density,  and  V  and  Vbaaerjum  are  the  volume  of  a  single  spore  and 
bacterium,  respectively  (eq  8,  Table  1).  A  single  spore  mass  (Table  1,  mg)  is  represented 

by  P spore ^ spore  =  m spore  ’  N spores  x  I03  is  the  number  of  spores  within  a  cubic  meter  (-  ,  )  ; 

m 

and  L  is  the  path-Length  ( m ).  The  effect  of  the  bacterium  and  spore  shape  is  taken  into 
account  when  computing  the  volumes,  V  and  Vhaclerjum  ,  in  eq  8. 

In  eq  1 1,  a  spherical  aggregate  with  a  lognormal  size  distribution, 
fn(r  :  p,(j),  was  assumed  as  given  in  eq  1,  with  a  total  number  of  aggregates,  N  ,  (in  a  unit 
volume)  and  a  total  volume,  V(p,cr),  of  the  aggregates,  shown  in  Figures  10  and  11. 

The  N  s  for  1  ACPLA  and  corresponding  CL  for  a  path-Length  of  100  m 
are  illustrated  in  Figures  8  to  15.  The  N  s  and  CL  for  a  family  of  lognormal  size 
distribution  fn(r  :  p,cr)  in  which  the  spore  and  bacterium  are  chosen  as  BG  are  presented  in 
Table  1.  In  Figures  8  to  1 1,  the  radii  range  in  fn(r  :  p,cr)  is  rmin  -»  0  and  rmax  =  25pm  .  In 
Figures  12  to  15,  rmm  — »  0  and  rmax  =  5 pm  were  selected  to  represent  a  scenario  in  which 
large  sized  particles  are  discarded  and  only  small  sized  particles  are  present. 

When  rmjn  =  0.5 pm  in  fn(r  :  p,cr),  the  radius  has  no  agent-containing 
particles  smaller  than  0.5pm .  The  CL  and  Nspores  computed  with  eq  11  are  higher  by  18% 

because  the  number  of  aerosol  particles  (N)  were  constructed  to  be  the  same  as  when 
rmm  —>  0 .  Thus,  the  size  distribution  fn(r  :  p,a)  will  contain  more  of  the  big  particles.  On 

the  other  hand,  if  rmin  — >  0  is  maintained  in  the  size  distribution  ( N)  in  eq  11  but  the  volume 
(V)  is  computed  only  for  rmin  >  0.5pm ,  the  value  for  N spores  and  CL  will  be  within  1%  of 
those  computed  with  rmin  -4  0  as  shown  in  Figures  16  to  23. 
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In  Figures  8  and  9,  N spores  and  CL  are  presented  for  Scenario  I  (dissemination  of 
bacteria).  In  Figures  10  and  1 1,  N  s  and  CL  are  shown  for  Scenario  II  (dissemination  of 

spores).  In  Figures  12  to  15,  the  situation  for  1  ACPLA  is  depicted  for  a  size  distribution 
with  only  small  particles  (r  <  5 jum).  Figures  12  and  13  show  the Nspores  and  CL  for 

Scenario  I  and  in  Figures  14  and  15,  Nspores  and  CL  are  presented  for  Scenario  II.  Figure  1 1 
provides  some  insight  into  the  crossover  of  the  curves  in  Figures  12  and  13. 


Figures  8  and  10  show  that  for  fn(r  :  /.i  =  l,cr  =  2)  (Fig.  1  b),  the  appropriate 
Af  „„  for  1  ACPLA  is  about  3  and  32  for  the  dissemination  of  bioaerosols  in  the  form  of 

spores 

bacteria  (Scenario  I)  and  spores  (Scenario  II),  respectively.  The Ns  s  for  the  dissemination 
of  only  spores  is  much  larger  (approximately  a  factor  of  10)  due  to  the  ratio 

[  —  -  -  ]  [  in  eq  11.  The  corresponding  CL  in  Figures  9  and  1 1  is  about  0. 1 1 

^ spore  ^ ’ bacteria 

and  1 .2  fug  /  m2 ,  respectively. 
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Figure  8.  Bacteria  or  Spore  N  es  Disseminated  as  Vegetative  Cells 
(calculated  with  eq  1 1  for  packing  efficiency  of  BG  bacteria). 
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Figure  9.  Corresponding  CL  to  Bacteria  or  Spore  Nspores  in  Figure  8 

(calculated  with  eq  1 1  for  packing  efficiency  of  BG  bacteria) 
(path-Length  =  1 00  m  N  =  ACPLA  =  1  l1). 


Figure  10.  BG  Bacteria  or  Spore  Count  (Table  1)  Disseminated  as  Only  Spores 
(ScenarioII)  (calculated  with  eq  1 1  for  packing  efficiency  of  BG  bacteria). 
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M  (M^n) 

Figure  1  1.  Same  as  Figure  9  but  Calculated  for^G  Spores. 


Table  2  displays  the  values  of  N spores  for  BG  dissemination  (Scenario  1  and 

II)  for  1000  ACPLA  and  the  corresponding  CL  values  for  a  path-Length  of  100  m  and  size 
distribution,  f„(r\ju  =  l,er  =  2)  (r  <  5  pm  and  r  <  25 /.on).  All  the  values  presented  in 

Table  2  can  be  easily  computed  with  eqs  1,  8,  10,  and  1 1,  values  from  Table  1  (first  2 
columns),  and  volumes  from  Figures  3  and  4.  The  estimated  CL  via  packing  efficiency  is 
smaller  than  the  estimated  CL  in  Section  4,  where  packing  efficiency  and  the  shape  of  the 
bacteria  were  not  taken  into  account.  The  CL  in  Table  2  is  also  smaller  than  the  estimated 
CL  given  via  the  “rule  of  thumb”  discussed  in  Section  6. 


Table  2.  Bacteria  or  Spore  Count  and  CL  for  100-m  path-Length  of  BG  Bioaerosols 
Disseminated  as  Bacteria  (Scenario  I)  and  Spores  Only  (Scenario  II). 


ACPLA  =  1000  (liter*') 

L  -  100  ( m ) 

BG  bacteria  and  spores 

f„{r  :/i  =  l,er  =  2) 

0  <  r  <  5  pm 

fn(r'P  =  h<7  =  2) 

0  <  r  <  25 jjm 

Scenario  I 

Scenario  II 

Scenario  I 

Scenario  II 

Number  of  bacteria  or 
spores 

N  N 

bacteria  ’  spores 

3022 

(~3  per  particle) 

32496 
(~30  per 
particle) 

4995 

(~5  per  particle) 

53712 
(-50  per 
particle) 

Concentration-path-Length 

CL^ f) 

m 

0.1092 

(—0. 1  pg/m2  per 
particle) 

1.174 
(~1  pg/m2 
per  particle) 

0.1805 

(-0.2  pg/m2  per 
particle) 

1.9405 

(~2  pg/nf  per 
particle) 
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Figure  12.  5G  Bacteria  or  Spore  Count  (Table  1)  Disseminated  as  Vegetative  Cells 
(Scenario  I)  with  Small  Particle  Size  Distribution. 
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Figure  13.  Corresponding  CL  to  Bacteria  or  Spore  Count  in  Figure  12  (Scenario  I). 
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Figure  14.  BG  N  t  (Table  1)  Disseminated  as  Spores  (Scenario  II)  with  Small  Particles 

Size  Distribution. 


Figure  15.  Corresponding  CL  to  N  es  in  Figure  14  (Scenario  II) 
(calculated  with  eq  1 1  for  packing  efficiency  of  BG  spores). 


9.  SUMMARY 

The  primary  objective  of  this  study  was  to  develop  relationships  and 
transformations  between  the  number  of  spores  (Nipores),  Agent-Containing-Particles  per 

Liter  of  Air  (ACPLA),  and  the  Concentration-path-Length  (CL)  for  aerosols  with  lognormal 
size  distribution.  For  the  estimation  of  N  s  and  CL,  2  methods  were  presented: 
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(I)  The  “rule  of  thumb”  (Section  6)  in  which  M  eK  was  taken  as  the  function  of  the 
radius  of  an  ACPLA  aerosol. 

(II)  Packing  efficiency  and  shape  superspheroid  for  bacteria  and  spores  (Section  8). 
For  method  II,  2  dissemination  scenarios  were  given: 

(a)  The  dissemination  and  aggregation  of  bacteria  to  form  a  lognormal  size 
distribution  fn(r:p,cr)  of  agent-containing  aggregates.  This  scenario  is 

applicable  to  non-sporulating  material  (e.g.,  Escherichia  coli). 

(b)  The  dissemination  and  aggregation  of  only  spores  to  form  the  lognormal 
size  distribution  fn(r:p,cr)  of  aggregates.  This  scenario  is  applicable 
for  sporulating  material  such  as  Anthrax. 

An  upper-limit  estimate  of  CL  (Section  4)  was  also  presented,  where  ACPLA 
particles  were  assumed  to  be  entirely  composed  of  spores  and  packing  efficiency  was 
ignored.  The  relationship  between  the  number  of  spores  and  CL  is  useful  for  computing 
fluorescent  cross-sections  from  bioaerosols.  For  example,  if  it  is  known  that  only  a  fraction 
of  the  spore  material  (e.g.,  Tryptophan)  is  the  source  of  the  sought  after  fluorescence,  the 
CL  estimate  can  be  adjusted  by  the  fraction  of  mass  due  to  the  material. 

The  estimate  of  N  s  and  CL  via  packing  efficiency  and  their  shape  was 
found  to  be  the  preferred  method.  This  study  has  provided  numerical  examples  for  Ns 

and  CL  for  1000  ACPLAs  and  a  path-Length  of  100  m  for  BG  bacteria  and  spores.  In  the 
examples,  fn(r  :  ju  =  l,cr  =  2)  was  selected  as  the  appropriate  aerosol  size  distribution,  as 

this  distribution  (Fig.  1  ab)  seems  to  be  a  viable  representation  of  BG  Bacillus  subtilis 
bioaerosols.  The  different  numerical  examples  are  presented  in  Table  3. 

Table  3.  BG  Bacteria  and  Spore  Estimations  by  3  Methods  of  Computation. 


ACPLA  -  1000  liter  1  with  aerosol  lognormal  size  distribution 
fn{r  :  /J  =  l,cr  =  2)  and  a  path-Length  of  100  m  of  BG  bioaerosols 

N 

spores 

(liter ) 

CL 

<*?) 

m~ 

Estimate  via  packing 
efficiency  and  effect 
of  bacteria  and  spores 
shape 
(Eqs.  8,11) 

Scenario  I 
(dissemination  ol 
bacteria) 

0  <  r  <  5  jdm 

3,000 

0.11 

0  <  r  <  25/am 

5,000 

0.18 

Scenario  II 
(dissemination  of 
spores) 

0  <  r  <  5  fim 

32,500 

1.17 

0  <  r  <  25 /.m 

54,000 

1.94 

Estimate  via  “rule  of 
thumb”  (Eq.  9) 

AWO  =  (  20255 

0  <  r  <  5/jm 

20,000 

0.71 

0  <  r  <  25  /am 

30,000 

1.07 

Upper  limit  estimate 
(Eq.  6) 

ACPLA  aerosols  are 
regarded  as  spores 

0  <  r <  5  pm 

2.6 

0  <  r  <  25  /am 

5.3 
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